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Abstract — It is well-known that adaptive homodyne estima- 
tion of continuously varying optical phase provides superior 
accuracy in the phase estimate as compared to adaptive or 
non-adaptive static estimation. However, most phase estimation 
schemes rely on precise knowledge of the underlying parameters 
of the system under measurement, and performance deterio- 
rates significantly with changes in these parameters; hence it 
is desired to develop robust estimation techniques immune to 
such uncertainties. In related works, we have already shown 
how adaptive homodyne estimation can be made robust to 
uncertainty in an underlying parameter of the phase varying as 
a simplistic Ornstein-Uhlenbeck stochastic noise process. Here, 
we demonstrate robust phase estimation for a more complicated 
resonant noise process using a guaranteed cost robust filter. 

I. INTRODUCTION 

Quantum phase estimation is the problem of estimating 
an unknown classical phase involved in the dynamics of a 
quantum system [1]. Precise phase estimation plays a key 
role in quantum computation [2], communication [3], [4] and 
metrology [5]. A fundamental bound on precision is imposed 
by quantum mechanics [6] and this limits gravitational 
wave detection [7] and can guarantee security in quantum 
cryptography [8]. 

Adaptive homodyne single-shot measurements of a fixed 
unknown phase can yield mean-square estimation errors 
below the standai'd quantum limit (SQL) [9], [10], [11], [12], 
[13], which is the minimum level of quantum noise obtained 
using standard measurements without real-time feedback. It 
is, however, practically more desirable to precisely track a 
continuously varying phase [14], [15], [16], [17], [18]. 

In Refs. [17], [18] the signal phase to be estimated is 
allowed to continuously evolve under the influence of an 
unmeasured classical stochastic Ornstein-Uhlenbeck (OU) 
noise process. However, since it is physically unreasonable 
to precisely know the underlying parameters of the noise 
process, the estimation process is significantly affected due 
to the uncertainty in these parameters. Hence, it is desired 
to make the phase estimation robust to uncertainties in these 
parameters. 

The authors have previously demonstrated in Ref. [19] 
that using a robust feedback filter designed based on a 
guaranteed cost robust filtering approach [20] can improve 
the estimation process as compared to an optimal filter In 
related works, the authors have also shown how continuous 
phase estimation using smoothing (rather than filtering alone) 
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can be made robust to uncertainties in the phase being 
measured and/or the measured output for a coherent beam 
of light [21] and a phase-squeezed beam of light [22], by 
employing robust fixed-interval smoothing theory from Ref. 
[23] for continuous uncertain systems admitting a certain 
integral quadratic constraint. 

These works, however, dealt with estimating the signal 
phase modulated by an Ornstein-Uhlenbeck noise process, 
which turns out to be a much simplified noise model as 
compared to the kind of noises that in practice corrupt the 
signal phase to be estimated. We, therefore, consider a more 
relevant and complicated second-order resonant noise pro- 
cess here and design a guaranteed cost robust feedback filter 
based on Ref. [20] for adaptive homodyne phase estimation 
of a coherent light beam phase-modulated by the same. 

II. RESONANT NOISE PROCESS 

The resonant noise process we consider here is the one 
typically generated by a piezo-electric transducer (PZT) 
driven by an input white noise as in Fig. [T] 

A. Transfer Function 

The simplified transfer function of a typical PZT is given 



by: 



G(s) 



(1) 



where k is the gain, is the damping factor, ujr is the 
resonant frequency (rad/s) and w is a zero-mean white 
Gaussian noise with unity amplitude. 

We use the following values for the parameters above: 
K = 1, C = 0.01 and ojr = 6.283 x 10^ rad/s, i.e. 1 kHz. 
The corresponding Bode plot of the transfer function dl} is 
given in Fig. |2] 

B. State-Space Realization 
From ([Til, we obtain: 

4>{t) + 2CuJr<j){t) + Ujl(j)(t) = Kv{t). 

Let xi :— (f) and X2 '■= (t>- Then, we have: 
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Fig. 1. Resonant noise 
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Fig. 2. Bode plot of the resonant noise process. 
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Let X : 

So, the state-space realization of the PZT output is: 



X = Ax + Gv, 



(5) 



where 



1 







where 

E[v{t)v{T)] = R(5(t- r), 
E[w{t)w{T)] = S(5(t-T), 
E[v{t)w{T)] ^ 0. 

Since v and w are unity amplitude white noise processes, 
both R and S are unity (scalars). 

The continuous-time algebraic Riccati equation to be 
solved to construct the steady-state Kalman filter for the 
system is then [24]: 

AP + PA^ + GRG^ PH^(JSJ^)^HP = 0. (10) 

The Kalman filter equation is [24]: 



i = (A - KH)x + KHx + KJio, 



(11) 



where K := PH^(JSJ^)-i is the Kalman gain. 

Using K = 1, C = 0.01 and Ur = 6.283 x 10^ rad/s, i.e. 1 
kHz as in section ITl-AI and |a| = 6 x 10^, we get: 

-10 
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and 



3.33785970 x 10"" 
8.02174133 X 10^^'' 



8.02174133 X 10" 
3.99751822 x 10"^ 



K 



4.80651797 x lO" 
1.15513075 X 10^ 



(12) 
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III. OPTIMAL KALMAN FILTER 

The standard adaptive homodyne phase estimation of a 
coherent state of light is optimal for given values of the pa- 
rameters when using a Kalman filter in the feedback loop to 
adapt the phase of the local oscillator Under a linearization 
approximation, the homodyne photocurrent from the adaptive 
phase estimation system is given by [17]: 



I{t)dt = 2\a\[(j){t) - 4>it)]dt + dW{t), 



(6) 



where |q;| is the amplitude of the coherent state with photon 
flux given by J\f :— jap, (f> is the intermediate phase 
estimate, and dW is Wiener noise arising from quantum 
vacuum fluctuations. 

The instantaneous estimate 6{t) is given by [17]: 



0{t) ^{t) 



m 

2\a\ 



= m 



1 

2W 



(7) 



where w := is a zero-mean Gaussian white noise with 
unity amplitude. 

Expressing (|7]i in terms of x defined in the previous 
section, we get the measurement model as: 



= Hx + Jw, 



(8) 



2|q 



where H := [ 1 ] and J 

Rewriting the equations for the system under considera- 
tion, we get: 



Process model: x = Ax + Gw, 
Measurement model: ^^ = Hx + Jw, 



(9) 



IV. ROBUST FILTER 

In this section, we make our filter robust to uncertainties 
in the parameters underlying the system matrix A using the 
guaranteed cost estimation robust filtering approach given in 
Ref. [20]. 

We introduce uncertainty in A as follows: 


where A := [ Si ] is an uncertain parameter satisfying 
||A|| < 1, i.e. Sl + Sj < 1, and < /ii < 1, < < 1 are 
parameters determining the levels of uncertainty. 

The process and measurement models from (|9) now take 
the form: 



Process model: x = (A + Di AEi)x + Gw, 
Measurement model; 9 = Hx + Jw, 



(14) 



Di 



El 



where 



As in Ref. [20], the Riccati equation to be solved in order 
to construct the guaranteed cost filter for the system is: 







eQEi^EiQ 



AQ + QA^ 

+ ^DiDi^ + GRG^ 



eQH^(eJSJ^ 
0. 



^HQ 



(15) 



The stabilising solution of this equation yields an upper 
bound Q for the robust filter error covariance. It is desired 
to obtain the optimum value of e (that we call egpt) at which 
the element in the first row and first column of the matrix 



Q (that we denote by Q+ := Q(l, 1)) is minimum. For 
example, for fii = 0.5, /i2 = and nominal values for other 
parameters, Fig.|3]shows the plot of versus e, where Copt 
is found to be 35. Thus, we get; 



^ 8.17703018 X 10"'" 4.09328251 x 10" 
The robust filter equation is: 



: (A + eQEi'^Ei - QH^(JSJ^)"iH)x 
QH^(JSJ^)"^Hx + QH^(JSJ^)"iJi( 



(16) 
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V. COMPARISON OF THE ROBUST FILTER WITH 
THE KALMAN FILTER 

A. Lyapunov Method 

We augment the system given by (fT4l l with the feedback 
filter ([TtI i and represent the augmented system by the state- 
space model: 

X = A X + B w, 



(18) 



where 

X := 

Thus, we have: 
A 



X 
X 



and 



V 

w 



A + DiAEi 
F L 



eQEi^Ei 



where F ;= QH^(JSJ^)-iH, L ;= A 
QH^(JSJ^) iH and 

G 

QH^(JSJ^) iJ 



For the continuous-time state-space model (fTsT i, the 
steady-state state covariance matrix P5 is obtained by solv- 
ing the Lyapunov equation: 



AP< 



— T T 

PsA +BB 



0, 



(19) 



where P5 is the symmetric matrix 



£;(xx^) 



Pi P2 

P2^ P3 



The state estimation error can be written as: 

e ;= X X = [1 — l]x, 

which is mean zero since all of the quantities determining e 
are mean zero. 

The error covariance matrix is then given as: 



S : = E{ee 
= [1 



= [1 - l]£(xx' 



1] 



Pi 



P2 
Pa 



Pi 



+ P3 



Since we are mainly interested in estimating xi = 0, the 
estimation error covariance of interest is cr^ = I](l, 1). 

B. Standard Quantum Limit 

The standard quantum limit is set by the minimum error in 
phase estimation that can be obtained using perfect hetero- 
dyne scheme. The SQL for our resonant noise process may 
be obtained using the same technique employed in Ref. [21] 
in the case of OU noise process. 

We use the fact that the heterodyne scheme of measure- 
ment is, in principle, equivalent to, and incurs the same 
noise penalty as, dual-homodyne scheme [17], such as the 
schematic depicted in Fig. |4l A signal at the input is phase- 
modulated using an electro-optic modulator (EOM) that is 
driven by the resonant noise source. The modulated signal is 
then split using a 50 — 50 beamsplitter into two arms each 
with a homodyne detector (HDl and HD2, respectively, with 
the local oscillator phase of HDl 7r/2 out of phase with that 
of HD2). The ratio of the output signals of the two arms 
goes to an arctan block, the output of which is fed into a 
low-pass filter (LPF). The phase estimation error would be 
minimum when this LPF is an optimal Kalman filter 

The output signals of the two arms are: 

Ii — —f=. (2|q;| sm(\) + ni + ^2) , 



V2 
1 

71 



{2\a\ COS0 + na — 714) , 



where ri i and 713 are measurement noises of the two homo- 
dyne detectors, respectively, and n2 and ri4 are the noises 
arising from the vacuum entering the empty port of the input 
beamsplitter corresponding to the two arms, respectively. All 
these noises are assumed to be zero-mean white Gaussian. 



Signal 




Fig. 3. versus e for fii = 0.5 and 112 = 0. 



Fig. 4. Block diagram of the dual-homodyne scheme for deducing the 
SQL for the resonant noise. 



The output of the arctan block is: 

2|ck| sin (/) + ni 



1} = arctan 



"-2 



2 a cos( 



(20) 



"-3 - "4 , 

Assuming the input noises are small, a Taylor series 
expansion upto first-order terms of the right-hand side yields: 

1 1 



—. — tTIi + — — -n2 



2\a 



2\a\ 



(21) 



Expressing this equation in terms of x, we get the mea- 
surement model as: 



?9 = Hx + Jz 



(22) 



1 
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where H = [ 1 ] and J = 

The error covariance matrix of the optimal steady-state 
Kalman filter for the process given by dD and the mea- 
surement given by ( |22] | may be obtained by solving an 
algebraic Riccati equation of the form (fTOt for P. The error 
covariance of interest (i.e. that in estimating xi — (j)) is then 
a2=P(l,l). 

C. Comparison of Estimation Errors 

The estimation errors may be calculated, as described in 
section I^A] for the robust filter, and likewise for the Kalman 
filter, as a function of 5i with 82 —Q for the nominal values 
of the parameters and chosen values for /ii with 112 = 0. 
These were used to generate plots of the errors versus Si to 
be able to compare the performance of the robust filter vis- 
a-vis the Kalman filter for the uncertain system with respect 
to the SQL. The SQL value is obtained by designing, as 
described in section IV-BI a different Kalman filter for each 
value of the uncertain parameter Figs. |5] |6]and|7]show these 
plots for III =0.2,0.5 and 0.8, respectively. 

Clearly, as Si deviates away from towards —1, the 
performance of the robust filter becomes superior than that 
of the Kalman filter in relation to the SQL for all levels of 

Ail- 
Similarly, the estimation errors may be calculated, and 
plots created, for the robust filter and the Kalman filter as 
a function of 62 with Si — for the nominal values of the 
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Fig. 6. Comparison of estimation eiTor covariance for the different 
filters witii fii = 0.5 and fi2 = 0. 
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Fig. 7. Comparison of estimation eiTor covariance for the different 
filters with /ii = 0.8 and 112 = 0. 



parameters and chosen values for /i2 with fii — 0. Figs. 
13 m and [To] show these plots for ^2 = 0.3,0.5 and 0.9, 
respectively. 

In this case too, the robust filter outperforms the Kalman 
filter as S2 tends towards —1. The SQL has not been shown 
in these plots since it is way above these errors. 
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Fig. 5. Comparison of estimation eiTor covaiiance a for the different Fig. 8. Comparison of estimation error covariance a for the different 



filters with fii = 0.2 and fi2 = 0. 



filters with fii = and fi2 = 0.3. 
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Fig. 9. Comparison of estimation eiTor covaiiance a'^ for the different 
filters witii /xi = and = 0.5. 




Fig. 10. Compaiison of estimation error covariance cr^ for the different 
filters with /xi = and ^2 = 0.9. 



VI. CONCLUSION 

This paper applies guaranteed cost robust filtering theory 
to the problem of robustly estimating the optical phase of a 
coherent state evolving continuously as a classical resonant 
noise process. We showed the behaviour of the robust filter 
as compared to the Kalman filter with uncertainties in the 
important parameters underlying the phase noise. The theory 
herein may be extended to include robust smoothing rather 
than filtering alone. Robustness to uncertainties in other 
parameters, such as the photon flux and noise power, may 
also be explored. Moreover, robustness for a squeezed state, 
instead of a coherent state, under the influence of such a 
resonant noise process would be interesting to investigate. 
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